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SUMMARY 



The present report deals with the flow p.ast an in- 
finite screen of thin airfoil { tv/o-dimens ional prohlem) . 
The vortex distrihution across the profile is established 
with appropriate expansion in series and the velocity 
distribution lift, .moipent,, and profile shape deduced. 
Inversely, the distribution is deduced from the vorticity. 

The method is the extension of the Birnbaum-Glauert 
method for the isolated wing. 

INTRODUCTION 



The present report describes a method for computing 
infinite screens of a.irfoils. 

This method is the natural but not immediate extension 
of the Birnbaum-Glauert method for isolated airfoils. By 
this is meant that the object of the present study involves 
screens of thin and slightly curved airfoils. 

For convenience the screen whose axis- (straigh:t line 
connecting the median points of the airfoil chord) is at 
right angle to the chord, is called straight screen ; that 
v^hose axis contains the profile chords, screen in tandem ; 
and that whose axis is oblique v/ithr the. chordJ, oblique 
screen . By the same definition the obliquity of the 
s creen is the angle of .the screen axis with the direction 
at right angle to the chord. 



♦"Sul caicolo. di schiere infinite di aili sottili. " Puhli- 
cazioni della R, Scuoia d • Ingegneria di. Pisa'j (sevefith- 
series) no,323, September 1937. 



2 



NACA Technical Me&iarattdtUD: No.. 968 



The term "chord" is employed iii :a Uroad^^sense and it 
signifies a straight line a little distance from the pro- 
file in reference to v/hich the ordinate© of the several 
points of the profile itself are measured, without impos- 
ing that it connect the extreme points or other conditdons 
of the kind. 



I. THE STRAIGHT SCREEN 



1. An infinite screen of eq.ual and equidistant vor- 
tices disposed on axis y, with circulation r and spac- 
ing h induces in a point P of the coordinate 2 = 
X + iy the velocity 

■°° i r ' ir rtz / , s 

w = - 2 ^ - — coth — (1) 

2iT(2-nih) ; 3h h 

corresponding to a complex potential 

0 = — In sinh 

2Tr • h 

(the axes are orientated as in fig. l). 
On a point of axis x we have: 

• = - in coth .m- (1«) 

2h . h . . . 

Denoting the components of the induced velocities by 
u and v., we put 

w = u ^ - . iv •' 
whence it is readily seen thsLt 

u = 0 V = |r coth 2£' (1«) 

2h h ■; 

Assume now that the vortex distribution follows a 
law Y s f(x) along the chord of the airfoils of the 
straight screen. 

Simply suppose that the wing chord is of length 21 
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and contains Taetween x = -I and x = I (fig. 2). Then 

the induced "veloislty of the whole screen on axis x is; 

I 

V ^ -i- Y dx' coth 3Li£z£lI (2) 
• 2h / h ■ . 



-I 



Then we put 



tanh ^ ^ i (3) 
h 



which readily yields 



coth HiLi) =iii£ 

\ h h / I - I' 



h 1 



dx = - : — ^ 



V 1- 



and hence 



V = ii- r Y 1 d|' (4)- 



where 



X = tanh 2L1. (5) 



We next put 



i ■ ■ 

- - cos 0 d I = X sine'd 6 (6) 



and find: 
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r p 1-7.^0 OS 9' COS 6* ' sin '6' .-y d 0' 
2tt / cos 9'- cosB l-X^cos^B' 

2tt / oos9'-cos 6 / . l-X^cos^ g' 

TT TT ■ 

" n Sine' 



_1_ 

2ti 



I / cos0'-cos9 - / 



■Vde'+^ • / sine'f 5^ Vae'l 

2 ,/ U-Acos.e« l+XcosG''^ J 



(7) 



It is interesting to note that the first of the tvro 
integrals of equation (?) is a function of 9, and con- 
seq.uently of x, while the second is independent or rep- 
resents a variation of the constant incidence. 

We then adop.t Otsor Y = f(6'). in Glauert's example, 
the series expansion 

V = V(ao cot -|- + S sin n 9') 

d 1 

where V indicates the asymptotic velocity along the 
positive direction of axis x, and , y^,... Yjj. . . , the 
single terms of the preceding series, and Vq . . . Vj^ . . . 
the corresponding values of v, so that v = 2 Vq. 

0 

In addition, '^ri ~ ''n ' + ''"n" where v^ ' is the 
part related to 9t v^^" the part not related to 9. 

V/e have* 

TT 0 1 TT . «« 

a.-, n sin9'cot-7r a^ p- 1+ cos 9' . ^ ^o 

27T / cos0»-cos9 Sir / cos9'-cos9 2 

■Jo Jq 



Vao^ r r^i+cos9' . i+cose' 



° " 4iT I / l-Xcos9''^° / l + Xcos9» 



0 



■J 

I 

cosh ~ tanh — = — ^^cosh — - l) 



♦For the calculus of the integrals in these formulas, see 

the appendix at the end of the report. 

**IoT. Vq , the first part Vq ' itself is independent of 
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whence y a, 



o ^ o * ' ^ T^"'-'= c 15*5 h ^ ( 8 ) 



Similarly, 



* t _ ^ r sin e« Bin n e"' , „, _ ' V an „„„ „ ^ 

V ' = -rr;- / — -— Q e ' = - — 5 — cos n 6 

^ 2tt / cos 9 ' - cos- e 2 

^ M _ IfM^ °r r ^sln ey sin n 6' , g ,_ /^ sin 9' sin n 6' , .,\ 
^ " 4tt Ij 1 cos 9' .y • 1 + X cos 9« J 



= coth Hi{tanh^ 111+ r-tanh HlV )= ^ 

4tt h 2h V -Shy J Ijan tanh^TTi ( 

2 2h ^ 



dd) 
is even) 



and hence 

y a, / V ni is odd 



^n = 



(- cos n 9 + ) (9) 

^ tanh^ ILL / n is oven 

2h 



And finally, ^sing 2* . to indicate that the sum is 
extended to odd values of n and S" to indicate the ex- 
tension to include even values: 



_ 00 

V* = Evji" = ^(a^j - E an cos n 9) 



V" =ILvn" = i 



tanh^ m 
2h 

2 a^ + I" an tanh^ |^ 

■ 1-tanh® ~ • 



- ■ -Sh 
or, by expansion in series: '* . 

• V" ^ | * 2 tanh^ |i . 

2. Mow the circulation corresponding to the vari- 

ous terms 7^ hence the complete circulation. T .= S 

can be computed. 



It i s t. 



Pn •= / Yn d X 
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where 



and therefore 



At k >^ sin 9 d 6 
■ ~ IT l-As cosse 



r = rY.pdx ^ ^ ^0 P 1 cos 6 cLe^Vao-hsinhlLi (lo) 

J ^ I l-X®cos3 e h ^ ' 



and ^ 



0 

TT 

p V a„ X h n sinn 0sin0 
r_ = / dx = — s / : i Q 

~ ^ ^ °- 0 (n is even) 



= V aj, h I 



tanh^ ^ (n is odd) 
3h 



(11) 



•whence 



00 

r = 2rn= Vh(^ao sinh^£^+ z'a^ tanh^ ^ ) 



y (12) 



= V h S ' (a„ + 2 a-,) tanh" Hi- 
1 " . 2h 

The velocity at infinity (negative for x = - », 
positive for x = +co, hut equal in hoth cases in magni- 
tude) is readily obtained from: 



CO 

2h 2 1'"^' ~~' 2h 



= ^ = 1 E'(aj^ + 2a) tanh° ^ (13) 



Froni the foregoing we can now draw some inferences 
ahout the differences between the direction of the tangent 
to the profile on the trailing edge and the direction of 
the velocity at <» , a difference which may "be termed angu - 
lar exaggeration of the blade screen. 

In effect, the trigonometric tangents of the two an- 
gles are expressed by viv)/Y. and Voa/V and their dif- 
ference 6 (which, knowing the assumed smallness of these 
angles, is dissipated in the difference of the angles) is 
expressed by 

-7 . 
where v(Tr) indicates the value of v for 6 = tt. 
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7 



Likewise, 8 can te considered as the sxun of 8^, 

each corresponding to T^; whence follows 

•nl 



^ I cosh Hi - sinh Lao ."Tr 
2 1 



h J 2 



odd 
even 



(16) 



3. If the flow velocity has a component -V sin a 
on axis y, the shape of the profile remaining th.e same, 
aQ may change in such a way that the Vq do not change. 
In other words, a change in a^ does not change the in- 
cidence of the profile. 

In the case where there is only and all the 

other terms are nullified, v ^puld also remain constant 
and the wing would be flat, with an incidence a ex- 
pressed hy 

a = cosh nl ' (17) 

2 h 

or 



a^ = ga ± (17') 

cosh 2l 
h 

and with the known formula (reference 1, p. 96). 

To = 2tt V I a ^Janh 2L1 (18) 

or, in comparison with the isolated wing, the circulation, 
and hence the lift, would he reduced in the ratio 



tanh 



Hi /ill 
h/ h 



In the general case equation (18) is still applicable, 
but the absolute incidence is expressed by 

' ■ ■ ■ ■ ■ " ■ - . \ ■ 

a „, ■ S< an tanh^ \ 
a = ^ cosh ILl + _i \ (19) 



tanh 

h 



8 
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We now have for lift P: 

P = p r y = 3Trp V^l ak (20) 
with k signifying the reduction ratio 

tanh 2i 

k = (21) 

Tf X 

h 

4. Proceeding to the determination of the shape of 
the wing, we have: 

= Z = afi. ^^gj^ El + 2" Vanh^ Hi - 2 ^ cos n 6 
dx ¥ 2 h 2 2h 2 

and so, as in the case of the isolated wing: 



an — "~ 



•1^ cos n e. d 0 (22) 



^ TT / dx 

'Jo 

It is understood that ~ should "be expressed in 

relation to 9 , 

Practically, given y ^and hence ^^as a function 

of X, we can pass from x: to 8 by means of the rela- 
tion 

, tanh IJ^ 

cos e = - ~ 

t anh "jj- 

Examule . - In the majority of cases the evaluation of 
the first three terms of the series will he sufficient to 
ohtain profiles of single and douhle camber.. 

For n = 0, it is simply 

dyp _ ao 

dx ~ 2 h 

which represents a straight line' of inclination ^ cosh 

For n = 1: 
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-. - -r^ .cos, e y ■ = - / cos 9 4f ^ ■ 

from which, after integration, fdllows 



, .■^ = a^ In cosh H£ + K 

^ ^ Trl ^ , TTl h 
2 tanh -35- 

Curve In cosh US. is shown in figure 3. 

h 

For n = 2: 

^ = ^ tanh« Hi - fli cos 2 9 
dx 2 2h 2 

which, integrated, gives: 



Vs = as 



tanh 5^ cosh ~ 

« X h 



ILl tanh^ nl ^ sinh^ Hi 
h h h 



The second term represents a constant inclination; 

the first, proportional to tanh — , has the shape of the 

curve in figure 3 and is repsonsihle for douhle camher of 
the profile. 

If the three terms are present, we have 
or else . .1 

Hi • 

.ml e~2h ag 



e ~h~ + ai- 



^ cosh # 2 cosh^ 

2h 2h 



for I = h, for instance, it is: 
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aestricted to the first term, it is: 

h , 

cosh ~- 
h 



The second reaches 0.041 ai and the third, -0.08 ag . 

5. We now make the comparison with the theory of suh- 
stitution vortices (reference 2). 

In the case of the flat wing (excluding s elf - induced 
velocity, disposing the vortices in the forward neutral 
point* and computing the induced velocity in the rear neu- 
tral point) it simply gives: 

p TT I 

v' = — coth - V a 

2h h 2tt I 

aixd on the other hand J 

r = - 2Trlv« 

Therefore: 

_r v_ ^^^^ y ^ 

2nl " i2h h ~ 2ttI ~ 



whence 

r = 2n V I a -~ tanh ^ 
ttI h 

which is the same formula previously worked out. Prom 
this it is concluded that the method of suhstitutidn vor- 
tices gives, for the flat wing, not merely approximate 
hut exact results. 

In the case of the caralDered wing (single camher) it 
is observed that the inclination of "both the leading and 
trailing edges ±ai/2. Now, in the case of the isolated 
wing, an equal inclination is obtained with an equal value 
of ai to w^hich corresponds a moment coefficient with re- 
spect to the forward neutral point, of ^ ai to which 



*The term "forward neutral point" and "rear neutral point 
for points situated at quart er--chord length from the tip 
was introduced by .K-assner. 
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11 



corresponds a dou'blet 



^J. = n...,.v,i' 



According to the 



general procedure of substitution vortices, it happens 
that we can compute the induced velocity in the rear neu- 
tral point of a screen of doublets (exclusive of that on 
the wing). But it is readily seen that the symmetry of 
the vortices with respect to the origin is such that the 
center of gravity of the circulation remains in the same 
origin for such a wing. In consequence, it is sufficient 
to consider the induced velocity of a screen of vortices 
disposed on axis y,' We will have: 



= ^ coth 



2h 



2h 



TTl 



Other than this, the form of the profile produces a 
change in the inc i a.eno e , or, what amounts to the same 
thing, a change in aero lift direction, given (for isolated 

wing) at As a result: 



= - 2Tr I r 



coth 



1- 



Hi 



2 h 



- TTl 



+ 2tt I V 



from which 



= V a^ h 



tanh 



Hi 
'2h 



^ 1 - A. tanh ILi 
•nl 2h 



A comparison with the exact formula 



aj^ h 



tanh 



TTl 

2h 



shows the approximation to he no longer satisfactory despite 
the fact that the ratio nl/h becomes of the order of 
unity or higher. 



It is interesting, on the other hand , to notice that 
the rigorous result can be obtained by assuming the induc- 
tion of screen vortices to be reduced to half. In effect 
then we have: 



r , = -TT I r 



1 



f.oth 



2 h 



1 

TTl 



+ 2tt I V 



from which the exact value for given above can be 

readily obtained. 
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6. A convenient and sufficiently approximate method for 
the practical prediction of the characteristics of a wing 
of given shape is that of stopping with the first three 
terms of the series., and deriving the coefficients of in- 
clination of the tangents to the profile in thre€ points - 
1, 0, I, "by the Birnhaum method. 

Denoting the inclination of the tangent in a generic 
point X, with pcj; and the inclinations . in points ABC 
(leading edge, trailing edge, and center of profile (fig. 
4) with a^, aj, aQ, vre find: 



°^A = 



aQ 



ao Trl 
— cosh — 



from which we obtain 



ao = 
As 



cosh ^ X cosh^~ 
h 2h 

- o^A 



ctA 
2 



+ + °-B 1 

2 cosh^Si 
2h 



and from eq.uation (l2): 



r = 7 h 



or, better: 
r = 2Tr 



tanh 



cosh 



ttI 



+ 2 tanh |1 
2h ■ 



JLi 1^ 2 * 2 



I 



(23) 



tanh^ 

2h 



^ (23») 
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In this very simple formula — = is the angle 



of lift for the isolated wing. The second term in paren- 
theses disappears when h/l — ^ w , while the reduction 

factor tanh ni. /ILL tends toward 1. 
h / 

Equation (23) or (23') can be used for all practical 

cases of the straight screen. 

It is noted that in place of equation (23) the fol- 
lowing procedure may he adopted: To compute r as the 
sum ©f Pq and each computed hy the method of suh- 

stitution vortices of which at n. 5, allocating to the 
profile incidence in the calculation of and 

a-Q — otp 

incidence — s ^ in the calculation of T, . It is 

2 

easy to see that the two methods lead to the same results, 

7. The prediction of the moment is more complicated. 
The moment M about origin 0 (center of profile) is ex- 
pressed "by 

M = pV / -Yxdx. 
- I 

but the expression of x obtained from equation (3) and 
equation (6) is too complicated for direct integration. 
It is more convenient to expand x in series of cosines. 

Posted, that is to say: 

CO 

X = E b^ cos n 6 
1 

we have : ^ 

= i r ^ cos n e de = - -2. r sin n 6 ^ __Sln i_ _ 



•J n 



Trn / IT l-A^oos^ 8 

o o 

which is easily obtained by integrating by parts; and, fi- 
nally, 

- tanh"" for n odd 

■h _ Trn 2h 
»■ ~ ^ 

0 f or n . even 
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00 

X = - -S^ S' i tanh^ ^ cob n 9 
and, consequently, v/ith the usual notation: 

TT 

M 



2 ■ r 

•^p A r a„ cot - i Z'^ tanh"' ^ cos m 6 

Jo 

TT i-Qo 

M_ rr-i^p a„ sinnG Z' tanh° 2l cos m8 

^0 

In particular, 

s 



sin 9d6 

1 - x^cos^e 

sin ed 6 

1 - x^cos^ e 



Mq = - p cosh 2i In cosh ^ (24) 



The point of application of the lift is ohtained "by- 
dividing Mq "by Pq = p 7® a^hsinh 2L1, whence: 



= - ^ coth . In cosh 211 (25) 
TT h h 



or else: 

tanh 2l1 



o . ^ fl + i tanh^ ILl + i tanh^ + ...) (25') 

° - \ 2 h S h • / 



2 Hi 



It is seen that Xq-> ~X for ^ -> 0 , whereas it tends 

tov/ard -1 for ■i--»co. 

h 

In other words, the aerodynamic center of the profile 
for the flat plate is shifted toward the leading edge so as 
to reduce the spacing (h) of the screen. 

As for the moments M^^ , we immediately have = 0 

for odd n, as iis shown hy the symmetry with respect to 0 
of the flow d is trihuti on , v;hile for even n the moment is 
other than zero, only the complete circulation and the lift 
heing zero. 



Mjj" is easily determined, when hearing in mind that 
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TT 



'Jo 



COS n 9 

i-A^cogs e 



d 6 for n < Q .is identical to that ob- 



tained for the same value of n with changed sign,. as 
distinguis'hed from that ohtained with the material appli- 
cation of the formula that gives the value of the said 



( 



integral I rr cosh — tanh^ — for even n 

^ h 2h . . • 



■) 



Therefore : 

3 



Mj, = - an?v S»'ltanh=i^ (Jm+n-i + ^m-n+i - Jm+n+l--rm-n-i ) 

2Tr 1 m 2h 

after posting (see appendix) 

TT 



Jn = 



cos n 



1- ?f cos^ 9 



d 9 



Thus we ohtain: 

, 3 



TT 



V3 a„ E' m tanh°^ 2ll . -3=- 



2h m 



tanh'^-^'^i itanh^-"" 2i 
2h 2h 



(26) 



where the + sign is valid for m ^ n - 1 and the - sign 
for m ^ n + 1. 

For n = 2, we have i 



M. = P V- a3 



( 



tanh^ ii- 



tanh 



a n I 
2h 



) 



in cosh ^+ 2(^: 



1+ tanh^ni^ 
2hy 



^^^Q V ag cosn -r— 

TT h 



2 In cosh 



TTl 



sinh» Hi 
h 



cosh« 



or, in more expressive form: 

2 r / 

Ma = — p V^ag tanh^ i + 2 tanh^ ( i+ 

T» 2h L 2h ^3 



tanh*2li tanh^3ll 
2h , 2h 



3 . 5 



5 • 7 



For 2-1 0, Mg -> 2-L- p ag, the value character- 
2h 4 ^ 



i6 
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istics.of the isolated wing, while for ^->«>, M2 — s> 0, 

Since a variation in incidence changes a_ only, 
the point Xq previously oDtained is such that with re- 
spect to it the moment does not change with -the incidence. 
This is-, moreover, the., aero dynam'ic center of " the • prof ile 
of the screen. 



II. THE SCHEEil,' IJST TANDEM 



8. The treatment is wholly similar to that of the 
straight screen. 

The induced velocity of a screen of vortices with 
equidistant h arranged on axis x (one of the vortices 
in the origin) is 

" i r 

-co 2tt(2 - n h) ~ 
and hence in a point on axis z: 



ir , TTZ 

cot — 

2h' h 



(27) 



- . V = -i- cot (27 «) 

2h h 

V/ith the same notation as "before, we post: 



tan ^ = I ; cotr^^ - V ^ ^ ' 



hence : 



V = 



TT 

JL r l+X^cos e cose' sin e ' y . a t roo'i 

2TT / cos 6' - cosa- ' 1+?nS.cos2 e' ^ ^ ^ 
d 0 



and, after simplv transformations: 



~ r sin e ' ^ ^0i_^ psin6'Cos8» 
Sir / cos e- ^ r cose ■ " ' 2ti ./ l+A^cos"* G*. 
o '"'0 
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wherefrom, after posting 



Y = V (ao cot I + S an sin n e) 



we olitain: 



Vjj' = - V cos ne 



an . .. I • 2 2h 



-V tan^ 5r for n even 



0 ■ for n odd 



and finally! 



v' = ( a« - Z a„ cos n 9 J (29) 



2 \ 0 ^ n 
as in the case of the straight screen 



= I bo (=03 ^ - .) . f « I-' .an= i^} 



or else 



V" = - i'^Can + 2 ao) tan^ 111 

2 2 n o 2h 

As regards the circulation, we have , correspondingly: 

• = V ao h sin - (30) 

r7- an h i^~^ tan" for ' n odd 

rn={ 2h (31) 

0 for n even 

and finally, for the flat wing: 

tan ^ 

■ s* 27t V I a • (S2) 

■ ■ • • TT I 

h • ■ • 
a formula similar to equation (18). 



18 
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The tandem arrangement amplifies, as is seen, the 
circulation considerably. 

Similarly, v;e can study th.e shape . of profiles corre- 
sponding to the various terms of the circulation, which 
is omitted here for the sake of "brevity. 

9. The comparison with the theory of suhstitution vor- 
tices remains to he made. For the straight wing, the 
vortices heing disposed in the forward neutral point , the 
induced velocity in the rear neutral point is: 

' t ^ * ^ V 

v» = — cot -——r - V a 

2h h 2ttI 

and 

r = - 2tt I v» 

whence : 

- JL. = il cot Hi - -I- - va 
2ttI 2h . h 2ttI 

from which we ohtain equation (32). 

In this case also, the results obtained are not 
merely approximate "but exact. And as for the straight 
screen, exact results are obtained for cambered v.'ing by 
disposing a screen of vortices in points in the center of 
the profile and reducing the induction by one-half. 

?or the sake of brevity, other similar considerations 
for the straight screen are omitted inasmuch as- regards 
the practical calculation of the lift in relation to in- • 
clination a^.aj , shape of profiles, and calculation 

•of the moment . 



III. TEB .OBLiqUE SCREEN- ■ 

10. Let 3 signify the obliquity of the screen, 
defined as the angle ..through which axis y* containing 
the vortices must rotate in order to be superimposed on 
axis y (fig. 6). 
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To pass from w referred to axes x' y', to w 
r ef e rred t o axe a-- X: y , — i-t i s su-f f 4 c-ien t • to po s-t in the 
complex potential 

, iP 
z * = z e 

;iPhus : ',. ■ 

w = - coth e^P . (33) 

2h h 

In a point on axis x, the value of w is pimply 

obtained by posing x at the pes ition 'of z. The value 
of w is complex, hence there is a real component, that 
is to say, along axis x, which is . negl acted , and an 
imaginary one, that is, along axis y, which alone is of 
interest, ■ 



We write : 



tanh I; A = tanh 33Ll§— 

h ' ■ h 



1- = - cos 6(6 complex) 
A 



(34) 



and merely follow the procedure step "by step, as for the 
straight screen. Thus 

■ = ? ao cot |- = V ao * ^ ' 

sin 9 

■Yn = V aji sin n 6 
with complex Rq . and ^q* Besides; 

ie^P r V ^ , V 4... TT(x-x»)e^^ 
w .= -. / Y dx* X coth 

2h / 

-I r ■ " ■ . 

whence 

- ^ _i__ P' 1-X cose cos 6' y . s in G ' " y ^ g, 

2n I cos 0'- cos 9 1- ?v'^cos^ 8' ' 
"--'o 

It follows that: 
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= ^ 5 Wo" - - (cosh ^^'^ - 1 J (35) 

and, on the whole i 

= IXaa cosh (36) 

Similarly: 

i T a / ° \ 

w = - cos n e + ■ ip \ ■ • (37) 

2 \ tanh'' nle__ ^ even 

\ .2h 7 

Proceeding with the calculation of T (generally com- 
plex) , we find : 

i 3 

= Vhe"^^ a^ sinh 2-^^ -• (38) 

n even 



""1 tanh'^' ILieif n odd 
L 2 h 



r„ = Vh e"^^ a„ J ifi (39) 



The total T must, of necessity, he real; in effect 
it contains an imaginary part which may represent a total 
source other than zero, which is naturally not admissihle, 
forcing the wing profile to he closed. Then we would have 

l{e"iP(ao sinh 2^ 4- S » a^ tanh^ Ili^ )) = 0 (40) 
L \ h " 2 h '''' 

where I. is the imaginary coefficient, 

11. The presence of a source distribution across the 
chord prompts us to attribute a certain thickness to the 
profile (E. Fistolesi - "Theory of Thin Airfoils" - in 
course of publication). We can, however, obtain the shape 
of the "dorsal spine" of the profile by the conventional 
method, namely, by ascertaining the vertical component of 
the induced velocity . t ' * 

▼ = f s{aQ cosh - ^ «^n cosne +2" a^tanh^lil^H } (41) 
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dy V 

where R Indicates the real part; and we put ~ T* 

For • yQ,' resulting from Vq, it yields the reoti- 

r : 

linear profile with incidence a = ^ R a- coth 2— 

. T 2 |_ ° h 

the yjj have curved profiles, from which the curve can 

also he computed; hut for simplicity this very complex 

development is foregone. 

If it is the only term in a^, we should have 



"o 

e^P 



s inh 



with Ap real, and then it follows that: 



2 0 " [_ 



a = -=-A- H e ^ coth 



(43) 



In general, the inclination of the profile for 6 = 
TT (trailing edge) can he calculated. It is: 

" i ^(^n cosh + 2 ' an - Z " an + 2 « an tanh° ^li^i 

V 2 \ . h 2 h 

„ r • cos s 

Por °>, we find: v^ = ^, or: 

2 h 



= i cos 3 {e"^^ (a-o sinh + 2t a,, tanh^ )} 

or else, since the part in hrackets is itself real; 

^ = I E (ao sinh Hijii ^ E ' a, tanh" ) 



whence the "angular exaggeration" 6 follows at: 
8 



4 , i■^ 

i'Tte . - . for n odd 



, r , , iB for n 

. i R (a^ e h~± Sa^ (l - tanh» Hi^i )) 



even 
(44) 
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12. A method for verifying that equation (40) is 
satisfactory is to put, as previously done in the case 
where a^ only is other than zero (equation (42)): 



.i3 



= -^o 



and in addition 



^n - 



s inh- 



TTle 



iP 



h 



> 



(45) 



tanh 
In this case: 



TT le 



i3 



(A-. real) 



2 h 



r= V h (A. + S« An) 



(46) 



In the case of the flat wing, the lift reduction fac- 
tor can he computed much more simply. In fact, posting 



r = 2TTlVak 

0 



where k is the reduction factor, leaves 



k ^ -fa 
a 21x1 



TT 



-E ( e^^ . coth 



h / 



(47) 



(48) 



and by simple calculation 

cosh^(^3^ cos 3^ - cos^ (^ILl sin 



rrl 



cos 



psinh (^ILlcosP^ cosh (^^^cosPj+sinPsin f^sinP^cos (^^^Lsinp) 

(49) 

or else: 



cosh 



/ 2Trl 



k = 



cos p^ - cos ^ sin p^ 



TTl 



cos P sinh 



( 



2lTl 



cos sin p sin 



( 



2tt1 



sin P 



) 



(49») 



h • y — ■ - — V b 

Formula (49) agrees with that hy Kumachi (reference 3), 
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accordiiig to a method previously IndioatecL hy Grammel 
(refference t) , with the difference, that the ,. Humachi- 
Grrammel procedure does not exactly f it 'the fliat wing; 
hut rather one with alight camher.* 

¥e can finally verify, for the case in point, the 
method of substitution vortices. 

The vertical induced velocity in the rear neutral 
point of the screen of vortices concentrated in the for- 
ward neutral point, the self-induced velocity, is 

r_R(eiP coth 

2h V h / 2ttI 

On the other hand, we have: 



from which equation (48) follows. Again the method of 
substitute vortices insures exact results. 

It will be noted that every change in incidence, 
without change in profile form, results in a change of 

; even the Aao corresponding to a change Aa must 
have the form (42), which must correspond to a zero lift 
source. However, in general, equation (47) holds true 
for r, or: 

r = 2-rT ; V-a k (50) 

the terms a^ , a^ , ag ...(and hence the shape of the pro- 
file) affect only the position of the zero lift curve 
vhose angle of lift a must be computed. 

The method, of substitute vortices can also be applied 
for obtaining Tj^ , if equation, (45) is taken for a^^ , 



*It should be remembered also that Grammel^s curved line 

is rather one-half the profile, the other half being 
formed by the symmetry of the first with respect to the 
chord; that it involves a flat rather than a curved pro- 
file of a certain thickness, which is more clearly shown 
in the treatment of the present study. 
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with, the single condition., that "the screen of vortices Tj^ 
are place.d-. in :the centers of :prof iles and the induced "• 
velocity is- halved'. 

'rfe have, in fact: 

i p 

r, = - TTir, I— z( e^^ coth H-i^ -i-")! + 2-iT I 7ai 

^ . . , ^ I 2h V .J8 h TTl/J 

where . aj,, the natural, i&cidence of the curveid profile, 
is given "by 

= 1 vx(0) -H ^x(^)., 1 H.(a,) 
^2 V 4 

from which follows: 

= M^il___vh 

R/^eiP eoth 
^ 2 h / 

or, with eq.uation (45) for ai, 

= 7 h Ai 

conf ormahle to equation (46) . 

■In general, the method of suhstitute vortices is ap- 
plicahle "by computing T as the sum of ■ and , 
giving the profile an incidence (fig. 4) for comput- 

ing Tq and incidence — z- for computing . It 

2 .» . 

affords 

... 

h 2h 

r . 2TrlV Jae — ^i— + 

H(e^3cosh-^) - H(e^^c.sh:ii|i^) 

(51) 
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13. The practical calculation of the aerodynamic 
characteristics of a given airfoil can he carried out hy 
stopping at the fifst three terms of the expansion of V, 

Equation (40) supplies i however: 

-ie . ^ Trle^P -iP ^ ^ TTle^^ \ 

I ( a„e sinh +■ a, e tanh ) = 0 

— V o h ^ 2 h / 

Moreover, with the notation used for the straight 
screen, we find: 



l1 

from which: 
a-B - o-A. - ECa^^ ) 

ac - "'^ ^ = R(a3) 

ac (l-R tanh^Hleii) tanh^ILlsli) = E(a^ coshUl^) 

J 

The prohlem generally remains indeterminate; to make 
it determinate two other conditions prevail (one, repre- 
sented hy equation (40)), which naturally must depend upon 
the profile thickness* 

Reserving the treatment of this point for the next 
paragraph, let us see what the results are from proceeding 
on the hasis of equation (45). We have-: 

•acfl-E tanh^ Hl^X 2ik±^(Y + E tanh^ 

A _ 2 h / 2 V 2 h y 

*o 



R 
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A ST ' I I .ll^JJl^^^Kr■yKP^«^^KaL.u^lMl■l■■^l> ' i i 



\ 2 h / 

and, substituting in equation (46): 

where k is expressed "by ei^uation (48), and 

2h 

= Hi ^ (54) 

H ^eiP cosh aLlsii") ' 
\ 2 h / 

the same k in which ^1 replaces length \, 

Equation (53/ coincides with (5l) when a^ = 0 and 
hence (a^ - a^,) + (a^ - Oq) = 0. Other forms into which 

equation (53) can he transformed are omitted for "brevity. 



14, It will "be noted that, with the effected positi.on, 
we have as a value of the source for x = 0, or for 0 = 

^, the imaginary coefficient of V \b.q cot S + ai sin , 
or : , 

I sinh 

\ ^ 

This source intensity defines an indiisat ion of the 

tangents to the profile 011 the "cop and tottom camher 
with respect to the tangents to the center line given by 

— , and, in consequence, an angle T hetv/een the two 
quoted tangents (fig., 7) given by* 

T = j(^e_i^ cosech Slj-li A^ l(^e^^ coth Hi^lij) (55) 



Such values of T . . are termed Tq* They vary, as is 
seen, w^ith the incidence, since A^ itself varies as the 
incidence varies. 



*Angl? T is positive when the nose is to the left. 
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The ensuing solution is a special case. In general. 



a- = (Aq + i Bo)e^^ cos^ch ""^^ 



i^ 



h 

aa = (Ajj + i Bn)e^^ coth" HJLgii 

2 h 

should replace equation (45). 

Vith these positions, equation (40) "becomes: 

Bq + S' Bj^ = 0 (56) 

while continuing the evaluation of equation (46), or: 

r = V h(AQ + Z' A^) (46) 

Stopping,, as usual, with the first three terms, and 
posting, for simplicity of notation. 



h 



e^P coth 2-Le = p + iq 



e ^ coth 2U:-S 

2 h 



> 



(57) 



e^^ coth"" JLlell = p^ + iq_ 
2 h 



(in which . k = — -fe — e ki— — \, we olitain in place of 
V -nip TT X -p^ / 

equation (52) the following: 

■ ' - = A ip 1 - B 3^q 1 ( 58 ) 

o-C ^-g ^ = -^sPs - Bs qs (59) 

.C (l - E tanh^ ^^i^)-^ m±i!^l+ R tanh** S^^^) = ^oP-Boq 

(60) 

Thes9 formulas are then supplemented by that ohtained 
"by angle T which, if it is considered that 
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cosech = coth - coth , is written as: 



h 2 h h 

From (56) finally. follows 



(61) 



+ = 0 



After -B.-jL is substituted for Bq , equations (58), 
(60), and (Sl) give three equations in A^ , 3^, and A^ 
which, resolved, give: 



Aq + Ai = (Aq + Ai)o + Bi M 
T = To + Bi N 



J 



(62) 



v/here (Aq + ^^.^ q a-^cL Tq indicate the values of Ag + A; 
and of T for Bj. = 3,,.= 0, anclj "besides: 



q,.. M + P - ^ 



1 



(63) 



As to T-, its value given in (55) can Tae written 



as : 



T, ={a,(l - E tanh^2liei!) 



TTleiP \ . ^B.-^ °^A 



(i + H tanh^ Jll)\ liZ 



+ (ttg a^) 



Pi 



2h/j p 
(64) 



For the profile of single •cani'ber, .so that ag = 0 



and hence a,Q = ■ — ^, it simply gives 



To = 2 ac ^i-^ + - 



(64») 



and finally: 
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r = (D- + V h ^ . AT . (-65) 

, ,. ._. ■» 

wrh.ere we p.ost . 

AT = T - To 

and (r)^ indicates the value of r for = 0, given 

in equation (53). 

This formula shows that T. depends not only on the 
form of the "mean line" of the profile, as for the straight 
screen and the screen in tandem, "but also on the airfoil 
thickness or, if preferred, to the curve of the top and 
hot torn camher . 

The prohlem is, as seen, rather complicated and the 
practical calculation of the otlique screen wich is the 
one that involves the propeller and the hlades of hydraulic 
machines, is necessarily la'borious. 



15. A few words on the moment calculation. The ele- 
mentary moment with respect to origin 0 is p T E(Yxdz) , 
and hence 

M = p7R/Yxdx 

To effect this integration, x is expanded in series 
hy posing co 

x= cos n9 

with _ 

rr TT 



i r X cos n 6d9 = ,-. r sin nS^ e'^P 9 dS 

W l-X^cos^e 

-'o 0 

whence f -ip 

. - ±£l^ tanh^ Trie 1 for n odd 

^n ■ -"i TT n 3 h 

0 for n even 

Procsediiig as for the straight screen, we find 

= -p hLvi H(a .-SiP cosh m cosh H^ei?.) 

°^ttV° h h'' 

If it is assumed that 

i P 

tto = A„ e^^ cosech ^^-^^ 



3^ • •N'AOA' -Teehaical' -MeitcJT^^^^^ 



the result is 
M = -p 



IT " \ h Ix / 



and now the point of application of the lift is character- 
ized hy : * ■ ■' ■ 



Similarly for M„: - 
= p E(a .-iPs^' tanh^ lOeii . i f tanh'^"''" Ilieii 



2 h ml 2 h 



±tanh»^-'"'lLl^j 
2 h i 



} 



where the + sign applies to m < n - 1, and the - sign 
to in ^ n + 1 . 



For 

■ an = A„ e^^ coth^ Ulsi?. 



the result is : 



} 
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APPENDIX 



The integrals used in the foreg^oing study are as 
follows: 



TT 



cos n9 ^- TTsinnBi 



... J Q _ _________ 

cos e - cos Oi " sin 6, 

The value of the' integral is that expressed "by the 

second term of the preceding formula even if 8 and 9i 
are complete TDecause the line along v/hich the integration 
is made, starting from 0 and terminating at tt, passes 
through the point 6^. 

3) Posting A = tanh a 

affords ^ 

T r cos n9 de ^ ^ ^n a 

l„ = / — = TT cosh a tanh — 

/ 1 - A cos 9 2 

TT 

T t r* cosn 6 d9 ^ / , a.^ 

I ' = / = TT cosh a( - tanh 2i ) 

jj, 1 + X cos e . ^ 2/ 

The ahove integrals reduced to the form 



f. 



cos n 6 d 8 tt p° 

l + p'^-Spcose .1-p' 



(See Laska, Collection of Formulas, p. 266), hy posing 

p = ^ -^}- - (p <C 1) = tanh % 

3) 

• TT r ^ 

« Q ^0 for n odd 

— °2i-o_e — d 9 = ^ 

1 - X® cos^ e ! TT cosh a tanh'^a for n even 

V 2 

Division of the integral hy two shows how easy it is 

to verify 

«3^n = i (In + In') 
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4) 



TT 

s in 9 s in n9 



o 



1 - cos® 6 



I- 



for n odd 

coth a tanh.^ S: for n odd 
2 



5) 

1 - X3 cos2 6 



TT 

COS 6 cos n a, ^ Q ^ 1 + 



for n even 



TT cosh^ a . T.n oo J, 

— ^ — — tanh" -K- for n odd 

l sinh a 2 



The above formulas for the integrals l^t ^'n» "^n* 
and their derivatives, are equally applicable in the complex 
field. 



Translation by J. Yanier, 
National Advisory Committee 
for Aeronautics. 



*Probably an error. 
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